Search for: All records

ABSTRACT High-velocity clouds (HVCs) may fuel future star formation in the Milky Way, but they must first survive their passage through the hot halo. While recent work has improved our understanding of the survival criterion for cloud–wind interactions, few observational comparisons exist that test this criterion. We therefore present an initial comparison of simulations with the Smith Cloud (SC; $d=$ 12.4 kpc, $$l, b = 40^{\circ }, -13^{\circ }$$) as mapped with the GALFA-HI (Galactic Arecibo L-Band Feed Array HI) survey. We use the SC’s observed properties to motivate simulations of comparable clouds in wind tunnel simulations with enzo-e, a magnetohydrodynamic code. For both observations and simulations, we generate moment maps, characterize turbulence through a projected first-order velocity structure function (VSF), and do the same for H i column density with a normalized autocovariance function. We explore how initial cloud conditions (such as radius, metallicity, thermal pressure, viewing angle, and distance) affect these statistics, demonstrating that the small-scale VSF is sensitive to cloud turbulence, while large scales depend on cloud bulk velocity and viewing angle. We find that some simulations reproduce key observational features (particularly the correlation between column density and velocity dispersion) but none match all observational probes at the same time (the large scales of the column density autocovariance is particularly challenging). We find that the simulated cloud (cloud C) showing growth via a turbulent radiative mixing layer (TRML) is the best match, implying the importance of TRML-mediated cooling for Milky Way HVCs. We conclude by suggesting improvements for simulations to better match observed HVCs. 

ABSTRACT The early growth of black holes (BHs) in atomic-cooling haloes is likely influenced by feedback on the surrounding gas. While the effects of radiative feedback are well-documented, mechanical feedback, particularly from active galactic nucleus (AGN) jets, has been comparatively less explored. Building on our previous work that examined the growth of a 100 $${\rm M_\odot }$$ BH in a constant density environment regulated by AGN jets, we expand the initial BH mass range from 1 to $$10^4\, {\rm M_\odot }$$ and adopt a more realistic density profile for atomic-cooling haloes. We reaffirm the validity of our analytic models for jet cocoon propagation and feedback regulation. We identify several critical radii – namely, the terminal radius of jet cocoon propagation, the isotropization radius of the jet cocoon, and the core radius of the atomic-cooling halo – that are crucial in determining BH growth given specific gas properties and jet feedback parameters. In a significant portion of the parameter space, our findings show that jet feedback substantially disrupts the halo’s core during the initial feedback episode, preventing BH growth beyond $$10^4 \, {\rm M_\odot }$$. Conversely, conditions characterized by low jet velocities and high gas densities enable sustained BH growth over extended periods. We provide a prediction for the BH mass growth as a function of time and feedback parameters. We found that, to form a supermassive BH ($$\gt 10^6 \, {\rm M_\odot }$$) within 1 Gyr entirely by accreting gas from an atomic-cooling halo, the jet energy feedback efficiency must be $$\lesssim 10^{-4} \dot{M}_{\rm BH} c^2$$ even if the seed BH mass is $$10^4 \, {\rm M_\odot }$$. 

Abstract The circumgalactic medium (CGM) of star-forming dwarf galaxies plays a key role in regulating the galactic baryonic cycle. We investigate how susceptible the CGM of dwarf satellite galaxies is to ram pressure stripping in Milky Way–like environments. In a suite of hydrodynamical wind tunnel simulations, we model an intermediate-mass dwarf satellite galaxy (M_*= 10^7.2M_⊙) with a multiphase interstellar medium (ISM;M_ISM= 10^7.9M_⊙) and CGM (M_CGM,vir= 10^8.5M_⊙) along two first-infall orbits to more than 500 Myr past pericenter of a Milky Way–like host. The spatial resolution is ∼79 pc in the star-forming ISM and 316−632 pc in the CGM. Our simulations show that the dwarf satellite CGM removal is fast and effective: more than 95% of the CGM mass is ram pressure stripped within a few hundred megayears, even under a weak ram pressure orbit where the ISM stripping is negligible. The conditions for CGM survival are consistent with the analytical halo gas stripping predictions in McCarthy et al. We also find that including the satellite CGM does not effectively shield its galaxy, and therefore the ISM stripping rate is unaffected. Our results imply that a dwarf galaxy CGM is unlikely to be detected in satellite galaxies; and that the star formation of gaseous dwarf satellites is likely devoid of replenishment from a CGM. 

ABSTRACT Arkenstone is a new scheme that allows multiphase, stellar feedback-driven winds to be included in coarse resolution cosmological simulations. The evolution of galactic winds and their subsequent impact on the circumgalactic medium are altered by exchanges of mass, energy, momentum, and metals between their component phases. These exchanges are governed by complex, small-scale physical processes that cannot be resolved in cosmological simulations. In this second presentation paper, we describe Arkenstone’s novel cloud particle approach for modelling unresolvable cool clouds entrained in hot, fast winds. This general framework allows models of the cloud–wind interaction, derived from state-of-the-art high-resolution simulations, to be applied in a large-scale context. In this work, we adopt a cloud evolution model that captures simultaneous cloud mass loss to and gain from the ambient hot phase via turbulent mixing and radiative cooling, respectively. We demonstrate the scheme using non-cosmological idealized simulations of a galaxy with a realistic circumgalactic medium component, using the arepo code. We show that the ability of a high-specific energy wind component to perform preventative feedback may be limited by heavy loading of cool clouds coupled into it. We demonstrate that the diverging evolution of clouds of initially differing masses leads to a complex velocity field for the cool phase and a cloud mass function that varies both spatially and temporally in a non-trivial manner. These latter two phenomena can manifest in the simulation because of our choice of a Lagrangian discretization of the cloud population, in contrast to other proposed schemes. 

Abstract The scaling of galaxy properties with halo mass suggests that feedback loops regulate star formation, but there is no consensus yet about how those feedback loops work. To help clarify discussions of galaxy-scale feedback, Paper I presented a very simple model for supernova feedback that it called the minimalist regulator model. This follow-up paper interprets that model and discusses its implications. The model itself is an accounting system that tracks all of the mass and energy associated with a halo’s circumgalactic baryons—the central galaxy’s atmosphere. Algebraic solutions for the equilibrium states of that model reveal that star formation in low-mass halos self-regulates primarily by expanding the atmospheres of those halos, ultimately resulting in stellar masses that are insensitive to the mass-loading properties of galactic winds. What matters most is the proportion of supernova energy that couples with circumgalactic gas. However, supernova feedback alone fails to expand galactic atmospheres in higher-mass halos. According to the minimalist regulator model, an atmospheric contraction crisis ensues, which may be what triggers strong black hole feedback. The model also predicts that circumgalactic medium properties emerging from cosmological simulations should depend largely on the specific energy of the outflows they produce, and we interpret the qualitative properties of several numerical simulations in light of that prediction. 

Abstract This paper presents a new framework for understanding the relationship between a galaxy and its circumgalactic medium (CGM). It focuses on howimbalancesbetween heating and cooling cause either expansion or contraction of the CGM. It does this by trackingallof the mass and energy associated with a halo’s baryons, including their gravitational potential energy, even if feedback has pushed some of those baryons beyond the halo’s virial radius. We show how a star-forming galaxy’s equilibrium state can be algebraically derived within the context of this framework, and we analyze how the equilibrium star formation rate depends on supernova feedback. We consider the consequences of varying the mass loading parameter ${\eta }_M\equiv {\dot{M}}_{\mathrm{wind}}/{\dot{M}}_{*}$relating a galaxy’s gas mass outflow rate ( ${\dot{M}}_{\mathrm{wind}}$) to its star formation rate ( ${\dot{M}}_{*}$) and obtain results that challenge common assumptions. In particular, we find that equilibrium star formation rates in low-mass galaxies are generally insensitive to mass loading, and when mass loading does matter, increasing it actually results inmorestar formation because more supernova energy is needed to resist atmospheric contraction. 

Abstract We present a suite of six high-resolution chemodynamical simulations of isolated galaxies, spanning observed disk-dominated environments on the star-forming main sequence, as well as quenched, bulge-dominated environments. We compare and contrast the physics driving star formation and stellar feedback among the galaxies, with a view to modeling these processes in cosmological simulations. We find that the mass loading of galactic outflows is coupled to the clustering of supernova explosions, which varies strongly with the rate of galactic rotation Ω =v_circ/Rvia the Toomre length, leading to smoother gas disks in the bulge-dominated galaxies. This sets an equation of state in the star-forming gas that also varies strongly with Ω, so that the bulge-dominated galaxies have higher midplane densities, lower velocity dispersions, and higher molecular gas fractions than their main-sequence counterparts. The star formation rate in five out of six galaxies is independent of Ω and is consistent with regulation by the midplane gas pressure alone. In the sixth galaxy, which has the most centrally concentrated bulge and thus the highest Ω, we reproduce dynamical suppression of the star formation efficiency in agreement with observations. This produces a transition away from pressure-regulated star formation. 

ABSTRACT In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. Active galactic nucleus jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work concludes from simulations of a halo with $$10^{14} \,\mathrm{ M}_\odot$$ that a successful jet model should have an energy flux comparable to the free-fall energy flux at the cooling radius and should inflate a sufficiently wide cocoon with a long enough cooling time. In this paper, we investigate three jet modes with constant fluxes satisfying the criteria, including high-temperature thermal jets, cosmic ray (CR)-dominant jets, and widely precessing kinetic jets in $$10^{12}-10^{15}\, {\rm M}_{\odot }$$ haloes using high-resolution, non-cosmological magnetohydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We find that scaling the jet energy according to the free-fall energy at the cooling radius can successfully suppress the cooling flows and quench galaxies without violating observational constraints. On the contrary, if we scale the energy flux based on the total cooling rate within the cooling radius, strong interstellar medium cooling dominates this scaling, resulting in a jet flux exceeding what is needed. Among the three jet types, the CR-dominant jet is most effective in suppressing cooling flows across all surveyed halo masses due to enhanced CR pressure support. We confirm that the criteria for a successful jet model work across a wider range, encompassing halo masses of $$10^{12}-10^{15} {\rm M_\odot }$$. 

Abstract Turbulent radiative mixing layers (TRMLs) form at the interface of cold, dense gas and hot, diffuse gas in motion with each other. TRMLs are ubiquitous in and around galaxies on a variety of scales, including galactic winds and the circumgalactic medium. They host the intermediate-temperature gases that are efficient in radiative cooling, thus playing a crucial role in controlling the cold gas supply, phase structure, and spectral features of galaxies. In this work, we develop an intuitive analytic 1.5-dimensional model for TRMLs that includes a simple parameterization of the effective turbulent conductivity and viscosity and a piecewise power-law cooling curve. Our analytic model reproduces the mass flux, total cooling, and phase structure of 3D simulations of TRMLs at a fraction of the computational cost. It also reveals essential insights into the physics of TRMLs, particularly the importance of the viscous dissipation of relative kinetic energy in balancing radiative cooling as the shear Mach number approaches unity. This dissipation takes place both in the intermediate-temperature phase, which reduces the enthalpy flux from the hot phase, and in the cold phase, which enhances radiative cooling. Additionally, our model provides a fast and easy way of computing the column density and surface brightness of TRMLs, which can be directly linked to observations. 

Abstract Traditional star formation subgrid models implemented in cosmological galaxy formation simulations, such as that of V. Springel & L. Hernquist (hereafter SH03), employ adjustable parameters to satisfy constraints measured in the local Universe. In recent years, however, theory and spatially resolved simulations of the turbulent, multiphase, star-forming interstellar medium (ISM) have begun to produce new first-principles models, which when fully developed can replace traditional subgrid prescriptions. This approach has advantages of being physically motivated and predictive rather than empirically tuned, and allowing for varying environmental conditions rather than being tied to local-Universe conditions. As a prototype of this new approach, by combining calibrations from the TIGRESS numerical framework with the pressure-regulated feedback-modulated (PRFM) theory, simple formulae can be obtained for both the gas depletion time and an effective equation of state. Considering galaxies in TNG50, we compare the “native” simulation outputs with postprocessed predictions from PRFM. At TNG50 resolution, the total midplane pressure is nearly equal to the total ISM weight, indicating that galaxies in TNG50 are close to satisfying vertical equilibrium. The measured gas scale height is also close to theoretical equilibrium predictions. The slopes of the effective equations of states are similar, but with effective velocity dispersion normalization from SH03 slightly larger than that from current TIGRESS simulations. Because of this and the decrease in PRFM feedback yield at high pressure, the PRFM model predicts shorter gas depletion times than the SH03 model at high densities and redshift. Our results represent a first step toward implementing new, numerically calibrated subgrid algorithms in cosmological galaxy formation simulations.